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Please add the following new claims: 
8. A method of determining a three dimensional shape, comprising: 
obtaining a first set of points; 

identifying local sets of points with associated vertices from the first set of points; 
forming surface equations representative of the local sets of points with the associated 
vertices; and 

ascertaining a plurality of curvature measures for the surface equations to identify a 



9. The method according to claim 8 wherein the surface equations are bivariate quadratic 
equations. 

10. The method according to claim 8 wherein the plurality of curvature measures is 
selected from the group consisting of: 

a first principal curvature; 
a second principal curvature; 
a mean curvature; 
a gaussian curvature; 
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a ratio of the minimum of the absolute value of the first and second principal 
curvatures to the maximum of the absolute values of the first and second principal curvatures; 

a difference between a first value of the first principal curvature and a second value of 
the second principal curvature; 

an average of the dot product between a first principal curvature vector of a surface 
vertex and a first principal curvature vector of all vertices on a 3D triangulated surface 
immediately connected to the surface vertex; 

an average of the dot product between a second principal curvature vector of a second 
surface vertex and a second principal curvature vector of all vertices on a 3D triangulated 
surface immediately connected to the second surface vertex; 

a difference between the average of the dot product between the second principal 
curvature vector of the second surface vertex and the second principal curvature vector of all 
vertices on the 3D triangulated surface immediately connected to the second surface vertex 
and the average of the dot product between the first principal curvature vector of the surface 
vertex and the first principal curvature vector of all vertices on the 3D triangulated surface 
immediately connected to the surface vertex; 

a radius or diameter of a circle inscribed into a normal triangle such that the three 
sides of the normal triangle are tangent to the inscribed circle; 

a radius or diameter of a circle circumscribed around a normal triangle such that each 
of the three vertices of the normal triangle intersects a perimeter of the circumscribed circle; 

a ratio of a radius or a diameter of a circle inscribed into a normal triangle such that 
the three sides of the normal triangle are tangent to the inscribed circle and a radius or a 
diameter of a circle circumscribed around a normal triangle such that each of the three 
vertices of the normal triangle intersects a perimeter of the circumscribed circle; 

an area of a normal triangle; 

a perimeter of a normal triangle; 

a ratio of an area over a perimeter squared of a normal triangle; and 

a dot product of a surface triangle normal vector and a normal triangle normal vector. 

11. A method of determining a shape, comprising: 
obtaining a first set of points; 

identifying a first predetermined number of local sets of points with associated 
vertices from the first set of points; 
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forming a first set of surface equations representative of the first local sets of points 
with the associated vertices; 

ascertaining a plurality of curvature measures for the first set of surface equations; 

and 

forming a shape equation with the plurality of curvature measures; 

12. The method according to claim 11 wherein the first set of surface equations are 
bivariate quadratic equations. 

13. The method according to claim 11 wherein the plurality of curvature measures is 
selected from the group consisting of: 

a first principal curvature; 
a second principal curvature; 
a mean curvature; 
a gaussian curvature; 

a ratio of the minimum of the absolute value of the first and second principal 
curvatures to the maximum of the absolute values of the first and second principal curvatures; 

a difference between a first value of the first principal curvature and a second value of 
the second principal curvature; 

an average of the dot product between a first principal curvature vector of a surface 
vertex and a first principal curvature vector of all vertices on a 3D triangulated surface 
immediately connected to the surface vertex; 

an average of the dot product between a second principal curvature vector of a second 
surface vertex and a second principal curvature vector of all vertices on a 3D triangulated 
surface immediately connected to the second surface vertex. 

a difference between the average of the dot product between the second principal 
curvature vector of the second surface vertex and the second principal curvature vector of all 
vertices on the 3D triangulated surface immediately connected to the second surface vertex 
and the average of the dot product between the first principal curvature vector of the surface 
vertex and the first principal curvature vector of all vertices on the 3D triangulated surface 
immediately connected to the surface vertex; 
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a radius or diameter of a circle inscribed into a normal triangle such that the three 
sides of the normal triangle are tangent to the inscribed circle; 

a radius or diameter of a circle circumscribed around a normal triangle such that each 
of the three vertices of the normal triangle intersects a perimeter of the circumscribed circle; 

a ratio of a radius or a diameter of a circle inscribed into a normal triangle such that 
the three sides of the normal triangle are tangent to the inscribed circle and a radius or a 
diameter of a circle circumscribed around a normal triangle such that each of the three 
vertices of the normal triangle intersects a perimeter of the circumscribed circle; 

an area of a normal triangle; 

a perimeter of a normal triangle; 

a ratio of an area over a perimeter squared of a normal triangle; and 

a dot product of a surface triangle normal vector and a normal triangle normal vector, 

14. The method according to claim 1 1 further including the steps of: 
obtaining a second set of points; 

identifying a second predetermined number of local sets of points with associated 
vertices from the second set of points; 

forming a second set of surface equations representative of the second local sets of 
points with the associated vertices; and 

determining values for the curvature measures for the local sets of points by the 
surface equations and applying the values to the shape equation to identify a shape. 

15. The method according to claim 14 wherein the second set of surface equations are 
bivariate quadratic equations. 
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REMARKS 

Claims 8-15 have been added. These claims are fully supported by the specification. 
It is respectfully submitted that no new matter has been added. Enclosed are the Replacement 
Claims Pages. 
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8. A method of determining a three dimensional shape, comprising: 
obtaining a first set of points; 

identifying local sets of points with associated vertices from the first set of points; 
forming surface equations representative of the local sets of points with the associated 
vertices; and 

ascertaining a plurality of curvature measures for the surface equations to identify a 

shape. 

9. The method according to claim 8 wherein the surface equations are bivariate quadratic 
equations. 

10. The method according to claim 8 wherein the plurality of curvature measures is 
selected from the group consisting of: 

a first principal curvature; 
a second principal curvature; 
a mean curvature; 
a gaussian curvature; 

a ratio of the minimum of the absolute value of the first and second principal 
curvatures to the maximum of the absolute values of the first and second principal curvatures; 

a difference between a first value of the first principal curvature and a second value of 
the second principal curvature; 

an average of the dot product between a first principal curvature vector of a surface 
vertex and a first principal curvature vector of all vertices on a 3D triangulated surface 
immediately connected to the surface vertex; 

an average of the dot product between a second principal curvature vector of a second 
surface vertex and a second principal curvature vector of all vertices on a 3D triangulated 
surface immediately connected to the second surface vertex; 
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a difference between the average of the dot product between the second principal 
curvature vector of the second surface vertex and the second principal curvature vector of all 
vertices on the 3D triangulated surface immediately connected to the second surface vertex 
and the average of the dot product between the first principal curvature vector of the surface 
vertex and the first principal curvature vector of all vertices on the 3D triangulated surface 
immediately connected to the surface vertex; 

a radius or diameter of a circle inscribed into a normal triangle such that the three 
sides of the normal triangle are tangent to the inscribed circle; 

a radius or diameter of a circle circumscribed around a normal triangle such that each 
of the three vertices of the normal triangle intersects a perimeter of the circumscribed circle; 

a ratio of a radius or a diameter of a circle inscribed into a normal triangle such that 
the three sides of the normal triangle are tangent to the inscribed circle and a radius or a 
diameter of a circle circumscribed around a normal triangle such that each of the three 
vertices of the normal triangle intersects a perimeter of the circumscribed circle; 

an area of a normal triangle; 

a perimeter of a normal triangle; 

a ratio of an area over a perimeter squared of a normal triangle; and 

a dot product of a surface triangle normal vector and a normal triangle normal vector. 

11. A method of determining a shape, comprising: 
obtaining a first set of points; 

identifying a first predetermined number of local sets of points with associated 
vertices from the first set of points; 

forming a first set of surface equations representative of the first local sets of points 
with the associated vertices; 

ascertaining a plurality of curvature measures for the first set of surface equations; 

and 

forming a shape equation with the plurality of curvature measures; 

12. The method according to claim 11 wherein the first set of surface equations are 
bivariate quadratic equations. 

13. The method according to claim 11 wherein the plurality of curvature measures is 
selected from the group consisting of: 
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a first principal curvature; 
a second principal curvature; 
a mean curvature; 
a gaussian curvature; 

a ratio of the minimum of the absolute value of the first and second principal 
curvatures to the maximum of the absolute values of the first and second principal curvatures; 

a difference between a first value of the first principal curvature and a second value of 
the second principal curvature; 

an average of the dot product between a first principal curvature vector of a surface 
vertex and a first principal curvature vector of all vertices on a 3D triangulated surface 
immediately connected to the surface vertex; 

an average of the dot product between a second principal curvature vector of a second 
surface vertex and a second principal curvature vector of all vertices on a 3D triangulated 
surface immediately connected to the second surface vertex. 

a difference between the average of the dot product between the second principal 
curvature vector of the second surface vertex and the second principal curvature vector of all 
vertices on the 3D triangulated surface immediately connected to the second surface vertex 
and the average of the dot product between the first principal curvature vector of the surface 
vertex and the first principal curvature vector of all vertices on the 3D triangulated surface 
immediately connected to the surface vertex; 

a radius or diameter of a circle inscribed into a normal triangle such that the three 
sides of the normal triangle are tangent to the inscribed circle; 

a radius or diameter of a circle circumscribed around a normal triangle such that each 
of the three vertices of the normal triangle intersects a perimeter of the circumscribed circle; 

a ratio of a radius or a diameter of a circle inscribed into a normal triangle such that 
the three sides of the normal triangle are tangent to the inscribed circle and a radius or a 
diameter of a circle circumscribed around a normal triangle such that each of the three 
vertices of the normal triangle intersects a perimeter of the circumscribed circle; 

an area of a normal triangle; 

a perimeter of a normal triangle; 

a ratio of an area over a perimeter squared of a normal triangle; and 

a dot product of a surface triangle normal vector and a normal triangle normal vector. 

14. The method according to claim 1 1 further including the steps of: 
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obtaining a second set of points; 

identifying a second predetermined number of local sets of points with associated 
vertices from the second set of points; 

forming a second set of surface equations representative of the second local sets of 
points with the associated vertices; and 

determining values for the curvature measures for the local sets of points by the 
surface equations and applying the values to the shape equation to identify a shape. 

15. The method according to claim 14 wherein the second set of surface equations are 
bivariate quadratic equations. 
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3-D SHAPE MEASUREMENTS USING STATISTICAL CURVATURE ANALYSIS 



DESCRIPTION 
BACKGROUND OF THE INVENTION 

Field of the Invention 

The present invention generally relates to a computer implemented 
method for applying differential computational geometry and statistics to 
detect shapes on three dimensional (3-D) computer models representative of a 
sequential series of two dimensional images of a structure. Twelve new 
5 curvature measures are also described. These measurements are applied to the 

surface of a 3-D computer model made up of triangles. Collections of 3-D 
images of known structures having aberrant surface curvature characteristics 
associated with dysfunction or defects are then used to train the software to 
recognize the particular curvature properties of an aberrant structure. Multiple 
1 0 linear regression is applied to the curvature measures of the 3-D test models 

and the results are optimized to remove multicolinearity and maintain a high 
coefficient of regression. 

Background Description 

An integral step in many medical diagnostic protocols is the detection 
15 of structural abnormalities associated with anatomical regions of interest. 

Accordingly, much of present day medical diagnosis involves the structural 
analysis of medical images either directly, through interpretation of film 
images, or indirectly, through the generation and interpretation of 3-D 
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computer models of an anatomical area of interest. Analysis of medical images 
can be difficult and time consuming, with the subjective view of the observer 
often clouding the interpretation of such images. A computerized method for 
providing objective measurements of the surface area and shape of an internal 
5 anatomical structure would therefore provide improved detection of 

abnormalities associated with structural defects. 

By way of example, cerebral aneurysms, also known as "berry 
aneurysms," are typically found as sacular outpoutchings of major intracranial 
arteries. They are known to enlarge progressively in many individuals and are 

1 0 thought to arise from arterial blood flow striking a weak area of the internal 

elastic lamina in the wall of an artery. Based on autopsy studies, these 
aneurysms occur in the general population at a rate anywhere between 1-8 %. 
For a given individual with an unruptured aneurysm, the annual risk of rupture 
is 1-2% while the lifetime risk of such a rupture is approximately 50%. At 

1 5 present time there do exist methods for elective treatment of such aneurysms in 

order to prevent rupture. Accordingly, effective diagnosis can potentially 
greatly reduce the morbidity and mortality associated with the occurrence of 
aneurysms. 

Traditional diagnostic protocols for detecting an aneurysms include 
20 cerebral angiography and magnetic resonance angiography (MRA). Cerebral 

angiography is an invasive procedure in which a catheter is advanced through 
the intracranial vessels and a contrast material is injected to enable visualization 
of the vessel walls. This procedure is costly, time consuming, and involves 
serious potential complications that can lead to aneurysm rupture, stroke, and 

2 5 death, although the risk of such complications is small (less than 1%). MRA, 

on the other hand, is less invasive, less time consuming, and presents less risk 
than cerebral angiography. 

Differential geometry has been used extensively as a means of analyzing 
a variety of geometric shapes and computational geometry is a means of 

3 0 rapidly determining surface characteristics using the mathematics of differential 

geometry. The techniques of computational geometry have been applied to the 
prediction of protein-protein interactions (Duncan et ah, "Shape Analysis of 
Molecular Surfaces," Biopolymers, 33: 231-238, 1993), analysis of dose-effect 
surfaces of combined agents (Lam, "The Combined Actions of Agents Using 
3 5 Differential Geometry of Dose-effect Surfaces," Bulletin of Mathematical 



Biology, 54:813-826, 1992), and the study of biological surface growth. 

Computational geometry has been applied in the medical field in the 
context of curvature analysis to diagnose lung cancer (Kawata et aL, 
"Computer Aided Diagnosis System for Lung Cancer Based on Helical CT 
Images " 1CIAP '97 Proceedings, 2: 420-427, 1997), cardiac aneurysms 
(Fantini et aL, "Quantitative Evaluation of Left Ventricular Shape in Anterior 
Aneurysms," Catheterization and Cardiovascular Diagnosis, 28: 295-300, 
1993), and assessment of the treatment of brain tumors (Dai et aL, 
"Intracranial Deformation Caused by Brain Tumors: Assessment of 3-D 
Surface by Magnetic Resonance Imaging," IEEE Transactions on Medical 
Imaging, 12, 4:693-702, 1993). All of these procedures have employed 
smoothing in order to avoid the large variations in curvature that are often seen 
when measuring the curvature of anatomic structures. The smoothing process, 
however, often distorts the measurement of curvatures. 

Recent developments have applied surface curvature analysis to extract 
feature lines on 3-D surfaces (Monga et aL, "From Partial Derivatives of 3-D 
Density Images to Ridge Lines/' Visualization in Biomedical Computing, 
VBC'92, Proceedings, pp. 118-129, 1992) but not as applied to 3-D surfaces 
of the cerebral vasculature. 

MRA studies of the brain commonly consist of 72-90 two dimensional 
(2-D) cross sectional images. Arteries in these images are highlighted using a 
3-D time-of-flight algorithm or its equivalent that is included as a matter of 
course with typical magnetic resonance scanner software. The images are 
displayed on sheets of 2-D radiological film for review on light boxes. Each 
image is read by a radiologist or surgeon who "creates" a three-dimensional 
interpretation of the vasculature through his or her knowledge of the anatomy 
of the blood vessels. This image may not be reduced to a three dimensional 
interpretation that can be reviewed by third parties. Due to the branching 
nature of the vasculature, this interpretation can be very difficult and time 
consuming. Branching vessels often appear as an aneurysm on a single image 
but analysis of surrounding images reveals normal vasculature. Similarly, an 
aneurysm might appear as a branching vessel unless it is carefully followed 
through a series of images to its termination. 

The use of computers to assist in detecting cerebral aneurysms prior to 
this invention has been limited. Attempts have been made to analyze image 
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data to detect cross-sectional images representative of aneurysms but these 
attempts have been futile, due to the difficulty in distinguishing an aneurysm 
from a normal branching point for an artery (Fessler et aL, "Object-Based 3-D 
Reconstruction of Arterial Trees from Magnetic Resonance Angiograms" 
5 IEEE Transactions on Medical Imaging, 10:1, 1991). Also, skeletal 

representations of the vasculature have been used to determine vessel branch 
points and topological features in vessels (Puig et ah, "An Interactive Cerebral 
Blood Vessel Exploration System" Visualization '97, Proceedings, pp.443- 
446, 1996). While these techniques are useful in some fields, such as virtual 

10 colonoscopy, the small vessel radii and large directional changes encountered 

in the topography of blood vessels are problematic to the application of 
skeletal representations for the identification of aneurysms. Smoothing and 
filtering of MRA images have been employed in an attempt to overcome this 
problem but the result is often a distorted view of the vasculature. 

1 5 Morphometric analysis has been proposed as a means of automatically 

detecting aneurysms but this procedure can only analyze small regions of the 
vasculature at a time (Matsutani, et al, Quantitative Vascular Shape Analysis 
for 3-D MR-Angiography Using Mathematical Morphology," Computer 
Vision, Virtual Reality and Robotics in Medicine, CVRMed '95, Proceedings, 

20 pp.449-454, 1995). 



SUMMARY OF THE INVENTION 



Based on the foregoing, a novel method for generating a more 
objective 3-D representation of anomalous regions on the surface of structures 

25 is disclosed. The present invention provides a method for objectively 

determining a 3-D analysis of a series of 2-D images. One particular 
embodiment of the invention is based on the simple observation that cerebral 
vessels are roughly cylindrical while aneurysms contained within such vessels 
are roughly spherical in shape. This particular embodiment exploits the 

3 0 spherical nature of aneurysms and computational differential geometry to 

identify and highlight aneurysms in unread 2-D MRA images. 

There are two classical mathematical entities used to analyze smooth 
surfaces which are known as the first and second fundamental forms. These 
measures are useful because they are intrinsic to the surface and therefore 

4 
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invariant to transform (rotation, translation, scaling). The mean and Gaussian 
curvatures are based upon the fundamental forms. 

In classical differential geometry, the principal curvatures k 1 and k 2? are 
averaged to produce the mean curvature (H) and multiplied to produce the 
5 Gaussian curvature (G). The mean and Gaussian curvatures can be used to 

classify the surface into particular types, In general, if H>0, a surface is 
concave and if H<0, the surface is convex. Beyond this, surface classification 
can be broken down into elliptic, hyperbolic, parabolic or planar according to 
the following rules. If G>0, the surface is elliptic, if GO, it is hyperbolic, if 

1 0 G=0 but H*0, the surface is parabolic and if G=0 and H=0, the surface is 

planar. It should be noted that cylindrical surfaces are parabolic and spherical 
surfaces are elliptic. Most algorithms that classify real surfaces according to 
curvature utilize an approach known as thresholding, wherein the principal 
curvatures are represented as a continuum and a particular threshold is used to 

1 5 separate the surface into different curvature types. 

One particular technique has been utilized previously to approximate 
the principal, mean and Gaussian curvatures at the vertices of a triangulated 
surface. This technique is based on the fact that a surface can be locally 
represented as a graph of a bivariate function and the curvature of this bivariate 

2 0 function can be directly measured (Hamann et al, "Curvature Approximation 

for Triangulated Surfaces," Computer SuppL 8; Geometric Modeling , 
Springer-Verlag Wein, New York, pp. 139-153, 1993). If the function 
provides a good approximation of the surface shape, then the curvature 
measured on the bivariate function will also be a good approximation of the 

2 5 underlying surface curvature. Since the interest in this process is the curvature 

at each vertex, a quadratic polynomial is fit to each platelet (region 
surrounding a vertex). 

The method of measuring curvature described previously (Hamman et 
aL) allows for a platelet to consist of only a single row of vertices around a 

30 central vertex (Figure 2). This single row of vertices is then used to determine 

the bivariate quadratic polynomial at each surface vertex. On "ideal surfaces" 
that are very smooth, this approximation works well. For surfaces that are not 
ideal, ones that have bumps and triangles of poor aspect ratio, this single row 
of surrounding vertices does not provide adequate data to determine a bivariate 
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quadratic polynomial that accurate approximates the local surface shape. 
Therefore, the technique described above has been modified to allow multiple 
rows of surface vertices to be added to the platelet and used in the 
determination of the bivariate quadratic polynomial The result is a much more 
5 accurate approximation of local surface shape and subsequently a more 

accurate approximation of surface curvature around each vertex. The 
determination of the number of additional rows of surface vertices to include is 
optimized by finding the number of rows of surface vertices that produces the 
largest coefficient of regression (described later). 

1 0 The classical curvature approximations provide a great deal of 

information about the surface shape and rate of change in surface shape. The 
classic measures are not, however, optimized to recognize and distinguish 
between cylinders and spheres. Surface size will also confound the application 
of classical curvature methods to cylinders and spheres. For example, a large 

1 5 spherical object may have the same Gaussian curvature as a small cylindrical 

object. Similarly, large and small objects of the same shape will have 
dramatically different curvatures, making it difficult to locate aneurysms, for 
example, with these methods independent of their radii. 

Twelve new curvature measures are disclosed that are useful to 

2 0 distinguish between spherical and cylindrical images presented in a 2-D format. 

Furthermore, these novel curvature measures are normalized to allow spheres 
to be recognized independent of their radii. Optimization of these measures 
through the application of multiple linear regression provides an optimum 
subset of curvatures that can be used to extrapolate structural anomalies from 

2 5 triangulated surfaces taken from 2-D image data sets. 

It is a further object of the invention to provide a method of evaluating 
determining three dimensional structures comprising the steps of: a) obtaining 
a computerized three dimensional representation of a structure or structures; 

b) identifying a first set of regions on the three dimensional structure or 

3 0 structures and assigning a numeric value to said structure or structures; 

c) identifying a second set of regions and assigning a numerical value to said 
regions; 

d) determining values for a plurality of curvature measures for each vertex on a 
surface of the structure or structures; e) performing multiple linear regression 

3 5 analysis on said values determined in said determining step to obtain a 
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coefficient of regression for all curvatures for all vertices; f) determining the 
variance inflation factor for each of said curvature measures; g) if all variance 
inflation factors are less than 10. go to step 1; h) if any variance inflation factor 
is greater than 10, sequentially reduce the subset of curvature measures used in 
5 multiple linear regression by 1 ; i) performing multiple linear regression on all 

combinations of curvature measures possible for each subset; j) selecting the 
subset of curvature measures that yields the largest coefficient of regression; k) 
performing multiple linear regression analysis on said values determined in d) 
to obtain a coefficient of regression for said curvature subset; and 1) inserting 
10 the partial coefficients of linear regression into the linear equation generated by 

multiple linear regression for said curvature subset. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The foregoing and other objects, aspects and advantages will be better 
understood from the following detailed description of a preferred embodiment 
15 of the invention with reference to the drawings, in which: 

Figure 1 is a flow chart representing a method in accordance with the 
present invention of applying novel and classical curvature measures and 
multiple linear regression to optimize structural differences in a series of 
sequential vertices. 

2 0 Figure 2 is a perspective view of the first principal curvature vector. 

Figure 3 is a perspective view indicating the normal vectors that 
comprise a normal triangle on the surface of a unit sphere. 

Figure 4 is a perspective view indicating surface normal vectors lying in 
a plane perpendicular to a cylinder. 

2 5 Figure 5 is a perspective view indicating the normal triangle 

characteristics for a sphere. 

Figure 6 is a perspective view showing the radius of an inscribed and a 
circumscribed circle in relation to the normal triangle. 

Figure 7 is a graphical representation of the unit normal vectors of a 

3 0 surface triangle and normal triangle. 

Figure 8 is a graphical representation of the correlation between the 
correlation coefficient and platelet radius. 

Figure 9 is a perspective view of a platelet constructed from vertices 
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DETAILED DESCRIPTION OF A PREFERRED 
EMBODIMENT OF THE INVENTION 

The present invention generally relates to a method and system as 
schematically represented in Figure 1, of applying differential computational 
geometry to the analysis of surface curvature. As shown in Figure 1, the 
method accepts as input a 3-D computer model of a structure that can be 
generated either automatically or semi-automatically from a collection of cross 
sectional images. 



Medical applications of the invention disclosed herein can be adapted to 
the study of much of the human body, as well as the bodies of other animals. 
Much of the human body is composed of organs that are roughly spherical or 
roughly cylindrical in shape. In many instances, derivation from one of these 
shapes is strongly correlated with disease. For example, the present techniques 
can be applied to the measurement of arterial stenosis, the detection of 
arteriovenous malformations, colonic polyps, and the detection of lung and 
liver cancers. 

The present invention can also be applied to the analysis of spherical 
and cylindrical machine parts that cannot otherwise be readily taken apart for 
inspection. In this particular embodiment, non-destructive testing of machinery 
via detecting anomalies in structure of internal parts of a machine. The 
invention is also applicable to the analysis of 3-D imagery generated from 
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weather patterns to detect funnel clouds or thunderstorms, for example. In 
another embodiment of the invention, the method can be applied to the 
technique of molecular modeling to detect abnormal protein-ligand interactions 
and in confocal microscopy detecting cancerous cells in a collection of images 
5 taken from a microscopic field or fields., 

Sequential cross sectional images of a selected anatomical structure are 
acquired with the use, for example, of a scanner such as a magnetic resonance 
imaging scanner (MRI). Any type of digital image scanner can be used in place 
of MRI images such as a helical computer tomography scanner (CT), 

1 0 ultrasound, or PET images. The 2-D images are arranged in a computer 

memory to create a 3-D data volume set. The image data to be analyzed can 
be generated and or stored in any of a variety of image formats. For example, 
the present invention is ideally suited for use with Picture Archiving and 
Communication system (PACS) format. The image data can be stored in the 

15 digital imaging and communications in medicine standard (DICOM), or as raw 

binary slices, or in a variety of volume formats. The image data can be stored 
in the computer memory in an internal data format which allows the image files 
to be saved as a single data volume instead of individual image files. The 
internal data format can also permit standard compression and uncompression 

2 0 techniques to be incorporated, reducing computer disk storage needs 



In one embodiment of the invention, a Web-based PACS software 
package acts as an automated filing system accepts and stores digital images 
created by traditional 2-D means. A second software component (IsoView) 
reads DICOM images stored in the PACS format and creates manifold (closely 

2 5 and singly connected )3-D triangulated surfaces from these images using a 

variant of the marching cubes algorithm. The surface images can be 
interactively displayed on a local workstation, stored as Virtual Reality 
Modeling Language (VRML) for later review, or stored as video recording or 
photographs for fiiture viewing. The 3-D images created in the IsoView 

3 0 format serve as the input for the curvature measurement method disclosed 

herein. 

The calculation of the bivariate quadratic polynomial at each surface 
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vertex is computed exactly as described by Hamman et al. The principal 
curvature measures, kl and k2, are then computed on these bivariate quadratic 
polynomials as described in the same. Hamman's technique determines the 
principal curvature magnitudes, kl and k2, but does not determine the principal 
curvature vectors, vl and v2. Therefore, equations and a mathematical method 
are described below to determine the principal curvature vectors vl and v2. 
The first step is to compute the eigenvectors (which were not described or 
computed by Hamman, et al.) of the same matrix for which the eigenvalues 
were determined, 

Equation 1 



C 2,0 C l,l 
C U C 0,2 



The eigenvectors are calculated by solving, 



k j " C 2 ? 0 



1,1 





V 




0" 








0 




j 







j = l,2 



Equation 2 



The solution to this system of equations is, 



- - e , K j - C 20 * 0, Equation 3 

J = l,2 



°u J 
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The 3-D principal curvature vectors vl and v2 are computed by normalizing 
the parametric principal curvature vectors and multiplying them times the 3D 
basis vectors determined in the method as described by Hamman, 




Equation 4 



The result are the two principal curvature vectors, vl and v2, that are 
used in newly disclosed curvature measures 3, 4, and 5 discussed below. 

As noted above, the classical curvature approximations provide a great 
deal of information about the surface shape and rate of change of surface 
shape, they are not optimized to recognize and distinguish cylinders from 
spheres. In addition, curvature magnitude also masks surface curvature 
analysis with the classical approximations. The present invention describes 
twelve novel scalar measures that are applicable to measuring curvatures. Two 
of these scalars combine the Gaussian and mean curvatures into a single value, 
three measure the change in principal curvature vector direction and seven 
make use of the surface normal vectors to predict surface shape. 

The calculation of these novel curvature measures are described as 
follows, where "||" represents the Euclidean norm: 

1. MIN(i|kl ||,||k2 ||)/MAX(l|Kl,|jk2 ||): A ratio of the two principal 
curvatures is calculated. Since both kl and k2 can assume positive 

2 0 (concave) or negative (convex) values, the unsigned magnitude of kl 

and k2 is determined. The smaller of the two numbers becomes the 
numerator and the larger becomes the denominator. This ration should 
be 0 for a cylinder, 1 for a sphere, and undefined for a plane. 

2. ||kl-k2|[ : The difference between the two principal curvatures is 
2 5 calculated. This value should be 0 for spheres and planes and 0 

otherwise. 
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Referring now to Figure 2, 

3. AVGOIv^v^l), j=l,n: The average of the dot product between the 
first principal curvature vector 200 at the center platelet vertex 201 and 
the principal curvature vectors 202 of those vertices immediately 

5 connected to this vertex is calculated. For convex surfaces, 200 will 

be parallel to the axia of a cylinder and will be randomly oriented on the 
surface of a sphere. Therefore this curvature measure should be small 
for cylinders and large for spheres. It also will be large in regions 
where the surface shape is changing dramatically such as at the branch 
1 0 point of vessels. 

4. AVG(lv 2 *V 2j :||) ? j=l,n: The average of the dot product between the 
second principal curvature vector at the center platelet vertex 201 and 
those vertices immediately connected to this vertex is calculated. For 
convex surfaces, v 2 will be perpendicular to the axia of a cylinder and 

1 5 will be randomly oriented on the surface of a sphere. Therefore this 

curvature measure should be large for spheres, small for large cylinders, 
and large for small cylinders. It should also be large at vessel branch 
points. 

5. AVGdv^v^ll), j=l,n) - AVG(lv 2 »v 2l :|j)J-l,n: The difference 
2 0 between the two curvatures in 3 and 4 above is taken. This measure 

should be small for small cylinders and large for large cylinders and 
spheres. 

Seven more curvatures disclosed herein use the surface normal vectors 
to predict surface shape. Referring now to Figure 3, all of these measures rely 
25 on one simple observation, i.e., for surface 300 as shown in Figure 3, if the 

normal vectors n 0 n, n 2 at the three corners of a triangle 301 which lies on 
surface 300 are drawn such that they originate from a single point, these 
normal vectors will form a triangle . For unit normal vectors n j? the vertices of 
this new triangle will lie on the surface of unit sphere 302. The circumference, 
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area and aspect ratio of this triangle 303 ( the "normal triangle") can provide 
considerable insight into the shape of the surface from which the triangle 
originated. Figure 5 illustrates the characteristic of the normal triangle for a 
cylindrical surface. A cylinder is defined by a circle extruded along a vector 
5 running perpendicular to the plane of the circle. Referring now to Figure 4, 

surface normals 401, 402, and 403 for a cylinder lie in planes perpendicular to 
the axis of the cylinder. If the unit normal vectors for three points on the 
surface of a cylinder are translated towards one another parallel to the axis of 
this cylinder, the three planes will eventually coincide and all three unit normal 

1 0 vectors will lie in this coincident plane. The "tips" of these unit normal vectors 

will form a circle 404 with a radius 1 unit larger than the radius of the cylinder 
405. If each unit normal vector is then translated in 2-D with the coincident 
plane, they can be positioned such that they originated from the same point 
but still existed entirely within the coincident plane. The tips of these unit 

1 5 normal vectors will therefore form triangle 406 that lies in the coincident plane. 

This normal triangle will typically have a large aspect ratio that will increase to 
infinity as the dot product of the normal vectors approaches 1.0. In other 
words, the more triangles on the surface of the cylinder, the larger the aspect 
ratio of the cylinder's normal triangles. Because the radius of the cylinder does 

2 0 not influence the aspect ratio of the normal triangle, this measure will be useful 

in detecting cylinders of any radius. 

Referring now to Figure 5, which illustrates the normal triangle 
characteristics for a sphere. All surface normals for a sphere are directed in a 
path from the center of the sphere through a surface vertex. If a plane is 

2 5 formed such that any two surface vertices and the sphere center are members 

of this plane, then the two unit normal vectors at the surface of the sphere will 
lie within this plane. The center of the sphere, the two surface points, and the 
tips of the two unit normal vectors will form two triangles 501 and 502, as 
shown in Figure 5. Equation 5, 

30 x^Cr+OVr 

defines, using the similar triangles theorem, the length of x m - one side of the 
normal triangle. If the same steps are carried out for the two remaining sides 
of the normal triangle, it becomes clear that each side is increased in proportion 
to (r+l)/r from the length of the surface triangle. If all three sides of a triangle 
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increase by the same proportion, the triangle formed (normal triangle) will have 
the same aspect ratio as the original surface triangle. Therefore the aspect 
ratio of the normal triangle depends on the aspect ratio of the surface triangle 
and is independent of the radius of the sphere, 
5 A better way to make use of the normal triangle curvature is to 

construct an equilateral triangle around the surface vertex of interest and find 
the normals at the vertices of this new triangle. Because the bivariate quadratic 
polynomial produces a good approximation of surface shape an equilateral 
triangle is created on this polynomial surface. Therefore the surface triangle 
1 0 area is uniform and does not influence normal triangle area. The normals at the 

vertices of this equilateral triangle are calculated using the quadratic 
polynomial coefficients shown below in Equation 6: 

*\ = °0,2 V + 0 1,1 U f * 



t* x I* 

il v 



1 5 The size, aspect ratio, and tilt of the normal triangle are used to predict 

the shape of the underlying triangles. For a spherical surface, the normal 
triangle aspect ratio will be the same as the surface triangle aspect ratio and 
will therefore also be an equilateral triangle, as shown in Figure 6. For a 
cylindrical surface, the normals will fan out in a straight line forming a triangle 

2 0 orthogonal to the axis of the cylinder. If the surface under examination is 

planar, all three normals will coincide and their tips will form a point. The size 
of the normal triangle also reflects the magnitude of curvature on the surface. 
Large normal triangles are seen in highly curved regions while small normal 
triangles are seen in relatively flat regions. Curvature measures that use these 

2 5 normal triangles are referred to as "normal triangle curvatures." 

The seven additional novel scalar measures are: 



6. Radius of an Inscribed Circle (inside the normal triangle) . (INSC) 
Referring now to Figure 6, there is only one way to fit a circle 602 
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inside the normal triangle 600 such that the three triangle sides are 
tangent to the circle and the circle center is at the intersection of the 
three angle bisectors. The three sides of the triangle have lengths L 0 , 
L„ and L 2 The radius of this circle 601 is calculated and should be large 
for a sphere, small for a cylinder, and 0 for a plane: 



S = -(L 



o 'Li 



L 2 ) 



Equation 7 



] 



(S-LXS-L.XS-L,) 



Equation 8 



10 



7. Radius of a Circumscribed Circle (CIRQ:Referring again to Figure 6, 
there is also only one way to fit a circle 603 through each of the normal 
triangle vertices such that the center is at the intersection of the 
perpendicular bisectors of the three sides of the normal triangle 600. 
The radius 604 of this circle 603 is calculated and should be very large 
for cylinders, 0 for planes, and small for spheres: 



15 



R 



circ 



— L 0 L 2 

4 R^c S 



Equation 9 



20 



8, Radius of an Inscribed Circle/Radius of a Circ umscribed Circle 
(RINCIR): This ratio measures the aspect ratio of the triangle. It 
should equal 0.5 for a perfect sphere and nearly 0 for a cylinder. This 
ratio is undefined for a plane. This ratio is relatively insensitive to 
curvature magnitude since this magnitude will appear in both th 
enumerator and denominator. It should therefore be large for spheres 
and small for cylinders independent of their respective radii. 

9. Area (NAREA): The area of the normal triangle is calculated by 
taking the cross product of two sides of this triangle and dividing by 

15 



2.0. It should be large for a sphere, small for a cylinder and 0 for a 
plane. 

1 0. Perimeter (NPER): The perimeter of a normal triangle should be 
large for a sphere, of intermediate length for a cylinder and 0 for a 
plane. 

1 1 . Area/Perimeter 2 (RNAP): Like measure 8, this measure is used to 
find the aspect ratio of the triangle. It is large for a sphere, small for a 
cylinder and undefined for a plane. Because the perimeter does not 
vary in magnitude as much as the circumscribed circle, this measure is 
more stable on cylindrical surfaces than measure 8. 

12. TILT : TILT is measured by first finding the unit normal vector of 
the surface triangle. This can be done by calculating the cross product 
of two of the three sides of the surface triangle and then normalizing 
(converting the vector to unit length) the normal vector of the surface 
triangle. The unit normal vector of the surface triangle is illustrated in 
Figure 7 as n s . Similarly, the unit normal vector, of the normal 
triangle is also determined. TILT is calculated by finding the dot 
product of these two vectors, n s and n nt This scalar value will be 1 for 
convex spherical surfaces, 0 for cylindrical surfaces and -1 for concave 
spherical surfaces. 

The above equations produce a group of statistical measures to 
determine curvature. A useful technique for combining multiple variables into 
a single equation is multiple linear regression. Multiple linear regression 
minimizes the sum of the squares of the residual between the regression 
equation and data to produce the linear equation that best fits the data. 
Furthermore, multiple linear regression uses the sum points over the entire 
surface. Significant variations in surface curvature ("Noise") seen in some 
curvature measures can be averaged out in this process. All that is required is 
for an underlying pattern to exist for the surface as a whole. Another attractive 



16 



WO 00/52629 



PCT7US00/05596 



feature of this technique is that the predictive ability of the equation improves 
as more 3-D surfaces are analyzed. In other words, the aneurysm predictions 
will become more accurate as more aneurysms are entered into the equation. 

Multiple linear regression is used to fit two or more variables to a set of 
5 data. The objective of the technique is to determine the coefficients to 

y=& + ftX, + &X 2 + ...)5„X m + e? Equation 10 

that produce the minimum sum of the squares of the residuals around the 
regression line as, 

S, =£(>', -ft -fa, -hh> A,)' Equation 11 

1=1 



10 In these equations, y is the dependent variable, x m are the independent variables 

and /3m are the partial coefficients of linear regression. Differentiating the 
above equation with respect to each of the regression coefficients results in the 
ma trix 
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Equation 12 



15 The x mj are the 16 curvature approximations consisting of the four classical 

curvature measures and the twelve novel measures reported herein, while y, are 
values assigned to the vertices of the 3-D computer models that 
4 teach"multiple linear regression what aneurysms look like. 

EXAMPLE 1 

20 In order to test the usefulness of the disclosed techniques in predicting 

aneurysms, MRA studies of 1 1 patients with normal arterial vasculature and 1 1 
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patients with diagnosed aneurysms were reconstructed. The 3-D models of all 
22 patients were created at a grayscale threshold of 325. The total number of 
vertices for all 3-D computer models is 290,802. The aneurysm patient group 
was contained nine females and two males while the normal group consisted of 
eight females and two males. The average age of the aneurysms group was 
53.2±17.7 years with a range from 22 to 78 years. The average age of the 
normal group was 44.6±15.3 years with a range of 24-73 years. 

After reconstruction, the surface vertices belonging to the aneurysm 
group were partitioned into three groups- aneurysm, transition region, or 
normal This was done interactively with IsoView. The entire dome of each 
aneurysm was partitioned and given a value of L0. The transition regions for 
each aneurysm were marked separately and given a value of 0.5. All other 
points on the model were assigned a value of 0.0. Ail vertices that make up 
the 3-D computer models of the normal patients were assigned a value of 0.0. 
These values constitute the y, dependent variables of the two previous 
equations. The transition regions lying between the aneurysm and normal 
vasculature is included because it has different surface curvature characteristics 
than the aneurysm dome and is often seen accompanying an aneurysm. The 
transition region is typically concave and is commonly referred to as the neck 
of the aneurysm. 

Having defined all variables, multiple linear regression of all data is then 
performed. The first step in this process is to compute all 16 scalar curvatures 
at each vertex of all 22 3-D computer models. The curvature approximations 
of equations 1-12 are measured for each model and stored in separate files for 
each patient and curvature type. A platelet radius determines the number of 
points to use in the least squares fit. The larger this radius, the more points are 
used. This modification solves the problem of directional bias that can occur 
when the image pixel size is much smaller than the distance between 
consecutive MRA images. Directional bias would cause the preferential 
selection of points in the direction normal to the image planes since triangles 
would tend to be "stretched" in this direction. No matter what radius is 
chosen, the first concentric row of points around the central vertex are always 
used. 

Because the regression of 16 variables over 290,802 vertices may be 
too large a problem to solve utilizing standard statistical software, software 

18 



was developed and optimized specifically to solve this large multiple linear 
regression problem. Such software was validated against previously solved 
problems found in Freund et al, Statistical Methods , Academic Press Ltd., San 
Diego, Ca, 1997 and Chapra et al. 9 Numerical Methods for Engineers with 
Personal Computer Applications , McGraw-Hill, Inc, New York, New York, 



A useful measure of the goodness fit of the regression line to the data is 
the correlation coefficient as in, 



and S t is the sum of the squares around the mean for the dependent variable y r 
The question then becomes what platelet radius to use in measuring the 
curvature on the 3-D computer models. The simplest way to answer this 
question is to try different radii to determine the radius that results in the 
largest correlation coefficient. 

The 16 scalar curvature measures for all 290,802 vertices were 
computed for platelet radii of 0.0 mm {i.e., one concentric row of points) to 
8.0 mm in increments of 1 .0 mm. Multiple linear regression was then 
performed on each group of curvatures and the regression coefficient was 
computed. Referring now to Figure 8, the relationship between correlation 
coefficient and platelet radius for the multiple linear regression of all 22 3-D 
computer models using all 1 6 scalar measures of curvature is shown. Clearly a 
platelet radius of 3.0 mm produces the best correlation coefficient: 0. 1 701 1 . 

Although this correlation coefficient does not indicate a perfect linear 
relationship, a test of the model was performed. The test for the model is 



1985, 




Equation 13 



where Sr is given by Equation 1 1 , 
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simply a test of the hypothesis that the entire set of coefficients associated with 
the 16 independent variables of equation (5.1) is 0. The alternative to this 
hypothesis is that one or more of the coefficients are not 0. The F statisticis 
defined by, 



2 

F(m,fl -m -1) = {n ~ mA)r Equation 14 

m{\~r ) 



is used to evaluate this hypothesis where r2 is the coefficient of determination 
(equal to the correlation coefficient to the power of 2), m is the number of 
1 0 curvature measures and n is the number of vertices analyzed. The null 

hypothesis is, 



H 0 :(pi,P2,...,p m) =0 - Equation 15 



The F statistic for our regression model using all vertices and all curvatures at 
a platelet radius of 3.0 mm is, 

15 F(16,290785)=541.55. 

The p=0.005 critical value for F(16 ~) is 2.14. Therefore, the null hypothesis 
that the regression model does not exist can be rejected and the partial 
coefficients p for the independent variables x are not all 0. The next step in 
confirming the model is a test of the individual partial regression coefficients. 

2 0 The test statistic for this hypothesis is also the F statistic but the equation to 

test a single partial coefficient f} } is, 

un iw ffi /e J> Equation 16 

F(ln-m-i) = — — — 
MSE 

where c^ is an element along the diagonal of the inverted multiple regression m 
x m matrix of equation (5.3) and MSE is the mean standard error defined by, 
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Equation 17 

MSE = ^ 



The null hypothesis is, 
Ho:j3=0. 

5 The p=0.005 critical value for F(l,«>) is 7.88. Table 1 lists the F statistic, the 

null hypothesis test and the partial (fi) and standardized Q3*) correlation 
coefficient for the multiple regression of all models and all curvatures for a 
platelet radius of 3.0 mm. The standardized partial correlation coefficients are 
defined by 

s 

10 p *=p %i Equation 18 

j ] s 



where s X] is the sample standard deviation of xi and sy is the standard deviation 
of y. These coefficients provide a means of comparing coefficients across 
variables by standardizing the means to 0 and variances to 1 . 

TABLE I. Hypothesis testing and partial regression coefficients for each 
1 5 curvature type. 
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Curvature 


F Statistic 


H 0 test 


A 


ft* 




N/A 


N/A 


0.01741 


N/A 


Kl 


-1.218 


cannot reject 
H 0 


-0.00459 


-0.08384 


K2 


-1.332 


cannot reject 
H 0 


0.01022 


0.22460 


G 


177.0 


reject H 0 


-0.00001 


-0.29593 


H 


-4.145 


cannot reject 

H 0 


-0.02413 


-0.10188 
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RK2K1 


0.002412 


cannot reject 
H 0 


0.00020 


0.000037 


DK1K2 


-0.05532 


cannot reject 
H 0 


-0.00083 


-0.03284 


DK1V 


-0.9754 


cannot reject 
H 0 


-0,9325 


0.10176 


DK2V 


-1.130 


cannot reject 
H 0 


-0.10054 


-0.14016 


DK1K2V 


-0.4760 


cannot reject 
H 0 


-0.06528 


-0.07455 


INSC 


0.4414 


cannot reject 
H 0 


-0.01578 


-0,08362 


CIRC 


68.08 


reject H 0 


0.01245 


0.03425 


RINCIR 


280.3 


reject H 0 


0.09167 


0.21795 


NAREA 


0.8571 


cannot reject 

TT 


0.00909 


0.00901 


NPER 


166.8 


reject H 0 


-0.01067 


-0.08601 


RNAP 


62.42 


reject H 0 


-0.82717 


-0.08710 


TILT 


1138.6 


reject H Q 


0.01030 


0.07563 



The above table discloses that the curvatures G, CIRC, RINCIR, 
NPER, RNAP, and TILT should contribute the most to the prediction of 

15 aneurysms. These six equations passed the null hypothesis. The preceding 

analyses demonstrate that a subset of curvature measures are useful in 
predicting aneurysms but does not provide any information as to how these 
variables are related to each other. Multiple linear regression attempts to 
determine the change in the dependent variable for a given change in an 

2 0 independent variable while holding all other independent variables constant. If 

two independent variables are colinear, it will be impossible to hold one 
constant while changing the other. The result of this violation is that one or 
more of the partial regression coefficients may attempt to describe a 
phenomena that is not exhibited by the data. This phenomenon is known as 
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multicolinearity, 

A useful measure of multicolinearity is the variance inflation factor 
(VIF), which indicates how much larger the variance of each independent 
variable is than it would be if it were not correlated with other independent 
5 variables. VIF is calculated by first performing a multiple linear correlation 

between each independent variable and all other independent variables. One of 
the independent variables x } is removed from the left hand side of equation 5 
and is used in place of y. The coefficient of determination (r j ) 2 is then 
determined for each independent variable and used to compute VIF: 

1 

VIF = — T 

1 ~ r j Equation 19 

10 



A VIF of 1 indicates that independent variable x, is not multicolinear 
with other independent variables. Anything larger than 1 indicates some 
degree of multicolinearity. Although there are no firm rules as to the 
magnitude that VIF can attain before multicolinearity is demonstrated, many 

1 5 statisticians have adopted a value of 1 0 as the VIF cutoff Thus, independent 

variables greater than 10 are considered multicolinear while those less than 10 
are considered independent. 

Table II lists the VIFs for the five classical curvature measures and the 
twelve novel measures of curvature described herein. Columns labeled subset 

2 0 size 1 6-5 represent the test of the correlation coefficient for a combination of 

16-5 combinations of curvature subsets. The curvatures were calculated on a 
bivariate quadratic patch created with a platelet radius of 3 .0 mm. Aside from 
CIRC and TILT, all curvatures are strongly multicolinear with each other when 
all 16 curvatures are analyzed. 

2 5 Having demonstrated a regression model, it must be determined which 

independent variables to remove in order to provide the most predictive model 
of curvature. This process is referred to as optimization. The endpoints of 
optimization are the maximum coefficient of regression with the minimum VIF, 
Because the coefficient of regression should decrease for each curvature 

3 0 measure removed, there is no single answer to this problem. In order to 
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determine the optimal subset of curvatures all possible subset combinations for 
a given subset size are examined through the process of an exhaustive search. 
The determination of the coefficient of regression for 2m subsets where m is 
the number of independent variables is called an exhaustive search. Multiple 
5 linear regression was performed on all 22 3-D computer models using all 

combinations of curvatures for a platelet radius of 3.0 mm. The results were 
sorted first by subset size (number of curvature measures used in the regression 
analysis) and then by the coefficient of regression. The curvatures that 
constitute the subset are those that produce the maximum coefficient of 

1 0 regression for a given subset size. Once the optimal subset is determined, the 

VIF for each curvature in the subset is computed. The coefficient of 
regression, F statistic and VIF for each subset is displayed in Table II. Empty 
entries indicate that the curvature measure was not used in the subset. 

As noted in Table II, there is no change in regression coefficient for 

15 subset sizes of eleven to sixteen subsets. Furthermore, different curvature 

measures can be included or excluded from an optimum subset, but RINCIR, 
NPER, and TILT are present as an optimal curvature measure for every subset 
from a size of three to sixteen curvature sets. 

In an other embodiment of the invention herein, optimization of subsets 

2 0 can be performed by the process of backward elimination, in which each 

independent variable is tested to determine the variable that can be eliminated 
while still maintaining the largest coefficient of regression. 
Table II indicates that curvatures H, DK1K2V, RINCIR, NPER> and TILT are 
ideal for identifying aneurysms when only five curvature measures are to be 

2 5 used. All five of these curvature measures were computed at each vertex of all 

22 3-D computer models using a platelet size of 3.0 mm. 
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Subset size 


16 


15 


14 


13 


12 


11 


Regression 
coefficient 


0.1701 


0.1701 


0.1701 


0.1701 


0.1701 


01701 


F statistic for model 


541.6 


577.7 


619.0 


666.6 


722.1 


787.7 




Kl 


190650 


212370 


192842 


217886 








K2 


1.98E+6 


1.40E+6 


1.40E+6 


39108 


205.3 


203.9 




G 


146.9 


146.6 


145.3 


1453 


146.6 


145.1 


r(VIF 


H 


34239 


35246 


36954 


7420 


25.03 


24.93 


RK2K1 


16.82 




16.49 


16.49 






o 
o 


DK1K2 


1.12E+6 


1.40E+6 


1.40E+6 








a 


DK1V 


2.10E+6 


3.36E+6 






1.310 


1.275 


§ 


DK2V 


2.40E+6 


3.36E+6 


2.10 


2.10 






re 


- DKIK2V 


2.I0E+6 


2.80E+6 


2.52 


2.52 


1.42 


1.40 


c 


INSC 


47.43 


46,50 






46.50 






CIRC 


5.16 


5.00 


5.14 


5.14 


5.00 


5.00 


> 


RINCIR 


50.70 


25.85 


48.24 


^ 48.24 


25.85 


25.08 


NAREA 


28.28 


26.15 


6.68 


6.68 


26.15 


6.01 




NPER 


13.32 


12.63 


9.36 


9.36 


12.63 


9.05 




RNAP 


36.45 


34.89 


23.59 


23.59 


34.89 


23.02 




TILT 


1.50 


1.50 


1.50 


1.50 


1.50 


1.50 



Subset size 


10 


9 


8 


7 


6 


5 


Regression 
coefficient 


0.1701 


0.1701 


0.1699 


0.1693 


0.1684 


0.1664 


F statistic for model 


866.2 


961.6 


1080 


1226 


1415 


1657 


Variance Inflation Factor (VIF) 


Ki 






110.7 


110.7 


110.6 




K2 


202.1 


194.9 










G 


143.7 


138.8 


110.2 


110.2 


110.2 




H 


24.91 


23.91 








2.12 


RK2K1 














DK1K2 














DK1V 


1.25 












DK2V 














DK1K2V 


1.38 


1.36 


1.36 


1.36 


1.36 


1.35 


INSC 














CIRC 


4.99 


4.99 


4.95 








RINCIR 


24.48 


24.16 


23.94 


18.87 


1.06 


1.05 


NAREA 














NPER 


5.52 


5.51 


3.46 


1.36 


1.34 


2.36 


KNAP 


20.26 


20.22 


19.83 


18.64 






TILT 


1.50 


1.47 


1 47 


1.46 


1.46 


1.07 
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Subset size 


4 


3 


2 


1 


Regression 
coefficient 


0.1655 


0.1624 


0.1538 


0.1203 


F statistic for model 


2046 


2625 


3522 


4273 




Kl 












K2 










r(VIF) 


G 










H 










RK2K1 






1.03 


1.000 




DK1K2 


1.35 








& 


DK1V 










B 


DK2V 












DK1K2V 












INSC 












CIRC 










/arian 


RINCIR 


1-05 


1.03 






TnFARBA 












NPER 


1.34 


l.Ol 








RNAP 












TILT 


1.03 


1.02 


1.03 





TABLE II (continued). 

At this point in the described method, the optimal curvature measures 
for the detection of aneurysms have been found through the optimization steps 
described above. The five curvature measures determined are not oniy optimal 
5 for the detection of aneurysms, but are also non-collinear. Therefore, each 

measure provides some useful information that can be used to recognize an 
aneurysm based on its shape. By performing multiple linear regression on these 
5 curvature measures for all vertices of all 22 3D computer models in this 
embodiment of the invention, it is possible to determine the coefficients, /3rn, of 
1 0 the equation 



y =0.017035 - 0.005998 H,+ 0.032173 DK1K2V 1 + 
0.04558 RINCIR, - 0.009688 NPER, + 0.012272 TILT r Equation 20 

The above equation is optimal for detecting aneurysms using the 
1 5 classical and new curvature measures described. The above equation can then 

be applied to any future 3-D computer models of aneurysms in the following 
way. A 3-D, triangulated computer model of the cerebral vasculature is 
created by any technique (e.g., marching cubes). A bivariate quadratic 
polynomial is fit to the collection of vertices, i.e., the platelet immediately 
20 surrounding a central vertex as shown in Figure 2. The subset of 5 curvature 

measures are calculated for each 3D surface vertex by determining their values 
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on the bivariate quadratic polynomial. Equation 20 is then used to determine y, 
for each vertex using the 5 individual curvature measures calculated for that 
vertex. A table of color values (a color table) is then created to represent the 
expected range of y 1 and the surface vertices or surface triangles are colored 

5 according to the color table. The color table can be interactively adjustable by 

the user or it can be fixed. Regardless of the technique used, the intention of 
the color table is to provide a simple means to color the 3-D computer model 
such that an observer's attention is directed to an area that has curvature 
properties consistent with an aneurysm. The 3-D computer models can be 

0 displayed in any number of different formats including 2D black-and-white 

images, 2D color images, stereographic images, movie format or in 3-D 
computer model format. The results of the above analysis can also be used to 
color or enhance the original 2D cross-sectional images to highlight aneurysms 
on the original 2-D MRA images, 

5 While the invention has been described in terms of a single preferred 

embodiment, those skilled in the art will recognize that the invention can be 
practiced with modification within the spirit and scope of the appended claims. 
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CLAIMS 

Having thus described our invention, what we claim as new and desire 
to secure by Letters Patent is as follows: 



1 LA method of evaluating determining three dimensional structures 

2 comprising the steps of: 

3 a) obtaining a computerized three dimensional representation of a 

4 structure or structures; 

5 b) identifying a first set of regions on the three dimensional structure or 

6 structures and assigning a numeric value to said structure or structures; 

7 c) identifying a second set of regions and assigning a numerical value to 

8 said regions; 

9 d) determining values for a plurality of curvature measures for each 

1 0 vertex on a surface of the structure or structures; 

11 e) performing multiple linear regression analysis on said values 

1 2 determined in said determining step to obtain a coefficient of regression for all 

1 3 curvatures for all vertices; 

14 f) determining the variance inflation factor for each of said curvature 

15 measures; 

16 g) if all variance inflation factors are less than 10, go to step 1; 

17 h) if any variance inflation factor is greater than 1 0, sequentially reduce 

1 8 the subset of curvature measures used in multiple linear regression by I ; 

19 i) performing multiple linear regression on all combinations of 
2 0 curvature measures possible for each subset; 

21 j) selecting the subset of curvature measures that yields the largest 

2 2 coefficient of regression; 

23 k) performing multiple linear regression analysis on said values 

2 4 determined in d) to obtain a coefficient of regression for said curvature subset; 

25 and 

26 1) inserting the partial coefficients of linear regression into the linear 
2 7 equation generated by multiple linear regression for said curvature subset. 
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1 2. The method of claim 1 wherein said plurality of curvature measures 

2 is selected from the group consisting of: 

3 the principal curvature kl ; 

4 the principal curvature k2; 

5 the mean curvature; 

6 the Gaussian curvature; 

7 a ratio of the minimum of the absolute value of the two principal 

8 curvatures to the maximum of the absolute values of the two principal 

9 curvatures; 

1 0 the difference between the first principal curvature value and the 

1 1 second principal curvature value; 

12 the average of the dot product between the first principal curvature 

1 3 vector of a surface vertex and the first principal curvature vector of all vertices 

14 on the 3-D triangulated surface immediately connected to this vertex; 

15 the average of the dot product between the second principal curvature 

1 6 vector of a surface vertex and the second principal curvature vector of all 

17 vertices on the 3D triangulated surface immediately connected to this vertex; 

1 8 the difference between the average of the dot product between the 

1 9 second principal curvature vector of a surface vertex and the second principal 
2 0 curvature vector of all vertices on the 3D triangulated surface immediately 

2 1 connected to this vertex and the average of the dot product between the first 

22 principal curvature vector of a surface vertex and the first principal curvature 
2 3 vector of all vertices on the 3D triangulated surface immediately connected to 
2 4 this vertex; 

2 5 the radius or diameter of a circle inscribed into the normal triangle such 

2 6 that the three sides of the normal triangle are tangent to the inscribed circle; 
2 7 the radius or diameter of a circle circumscribed around the normal 

2 8 triangle such that each of the three vertices of the normal triangle intersects the 

2 9 perimeter of the circumscribed circle; 

3 0 the ratio of the radius or diameter of a circle inscribed into a normal 

3 1 triangle such that the three sides of the normal triangle are tangent to the 

32 inscribed circle and the radius or diameter of a circle circumscribed around a 
3 3 normal triangle such that each of the three vertices of the normal triangle 

3 4 intersects the perimeter of the circumscribed circle; 
3 5 area of the normal triangle; 
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3 6 perimeter of the normal triangle 

3 7 the ratio of the area over perimeter squared of the normal triangle; and 

3 8 the dot product of the surface triangle normal vector and the normal 

3 9 triangle normal vector. 

1 3. The method of claim 1 wherein said three dimensional representation 

2 is obtained from an image of said structure. 

1 4. The method of claim 2 wherein said image is constructed from 

2 multiple two-dimensional scans. 

1 5, The method of claim 1 wherein said three dimensional structure is a 

2 tumor, aneurysm, or polyp, and said predetermined shape is spherical. 

3 6. A method of evaluating three dimensional renderings of blood 

4 vessels for detecting the presence of aneurysms comprising the steps of: 

5 a)obtaining a computerized three dimensional representation of said 

6 blood vessels;; 

7 b) identifying a first set of regions on the three dimensional structure or 

8 structures and assigning a numeric value to said blood vessels; 

9 c) identifying a second set of regions and assigning a numerical value to 

10 said regions; 

11 d) determining values for a plurality of curvature measures for each 

12 vertex on a surface of said blood vessels; 

13 e) performing multiple linear regression analysis on said values 

1 4 determined in said determining step to obtain a coefficient of regression for all 

1 5 curvatures for all vertices; 

16 f) determining the variance inflation factor for each of said curvature 

17 measures; 

18 g) if all variance inflation factors are less than 1 0, go to step 1; 

19 h) if any variance inflation factor is greater than 10, sequentially reduce 
2 0 the subset of curvature measures used in multiple linear regression by 1 

21 i) performing multiple linear regression on all combinations of 

22 curvature measures possible for each subset; 

2 3 j)selecting the subset of curvature measures that yields the largest 
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2 4 coefficient of regression; 

25 k) performing multiple linear regression analysis on said values 

2 6 determined in d) to obtain a coefficient of regression for said curvature 

2 7 subset; and 

28 1) inserting the partial coefficients of linear regression into the linear 

2 9 equation generated by multiple linear regression for said curvature subset. 

1 7. The method of claim 1 or claim 6 further comprising the steps of; 

2 determining a single scalar value for each vertex on a structure or 

3 structures from the linear equation determined; 

4 creating a color table of all scalar values for each said vertex; 

5 assigning a color to each vertex; 

6 rendering the 3-D structure or structures with said color values 

7 displayed on the surface of the structure or structures. 
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FIG.1 
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DETERMINE VALUES FOR A PLURALITY OF CURVATURE MEASURES 
FOR EACH VERTEX ON A SURFACE OF THE STRUCTURE 



PERFORM MULTIPLE LINEAR REGRESSION ANALYSIS TO OBTAIN 
A COEFFICIENT OF REGRESSION FOR ALL CURVATURES FOR ALL VERTICES 
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CURVATURE MEASURES 
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REDUCE THE SUBSET OF CURVATURE BY 1 



PERFORM MULTIPLE LINEAR REGRESSION ON All COMBINATIONS 
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TO OBTAIN A COEFFICIENT OF REGESSION 
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